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Navier- Stokes Computations on Full Wing Body 
Configuration with Oscillating Control Surfaces 

coefficients with and without contro s “« a e J^ contro , surfaces are modeled by usmg moving gr.ds^ 
experiment. For unsteady cases, oscillating traih g * dilatory motions of the control surfaces on 

Response characteristics between symmetric and hif Lr lift than the steady case with no 

the left and right wings are studied. The ant **^ e h |ift than the antisymmetric case. The detaile 
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Introduction 

conec. {' SSL of <1,0 aircraft s,a- 

rb” ,hc inves ' 

BSSS^SSS 

Furthermore, aircraft are often J uns t ea dy aero- 

lation because of the flow unsteadin - h vily on 

dynamic environment, many modern jT'? er Sn 
active controls for safe an ^ b ® cn studied as a design 
An arrow-wing confl f“'.f r , . Because of the highly 

concept for supersonic u P ic flutter is a design 

swept thin wing, it is know ‘ opment of an analytical 
problem on this configura it ■ aeroelastic performance of 

recently developed 1 ci • S simultaneously integral- 
computing aeroelastic response > modal struc- 

SSSsSS 

free license w esenase ali h ‘ “ ,^. rvcd by ,hc 

for Governmental purposes. All other rig 

tfc; S^2sS!Sd£?SS^S^. 


wings undergoing un^e^d^motions^ '^e code ^ osci u at ing 

to simulate unsteadyflows ov^jdw g ^ ^ 

- Hand, a fu.l-span 

Jing-body configuration ^.“J^iS^avier-Stokes 

s Ju“o^ 

Computed pressures and g ■ eKoeriment. 1 Compar- 

ten compared w.th '* and an- 

JiT— "rface mo.ions on ,h« right and left wines 

is also presented. 

Numerical Method 

The nondime nsionalized Th"e vSoshj 

lisISSSSH: 

mmrnm 

mmm 

wmm 

single Cray-C90 processor (175 MFLOFts ana h- 
single Cray Y-MP processor). 

Model Geometry and Grid 
T « II m toD0 loav grid is used for a wing-body configu- 
,aS ”lh aZ,S,<ace Th« topotopy is chosen m order 
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f K ' VVlnd ‘ Unnel model 8«*metry of an arrow- wing configuration (all dimensions in centimeters). 



Fig. 2 Overview of the surface grid. 

DD e N S rA a rl ign f t grid lmeS *° ' he COntro1 surface ' T he ICEM 
^ Fmm ,h Ur / SyStCn r Was used to ^nerate the surface 
grid. From the surface grid, the volume grid was generated 

hasTwo fl ins T^Ncode. 1 '' A| th°ugh the experimental model 1 
as two flaps both at the leading and trailing edges onlv the 

S ,n,il,ng edge “ COnsidered in 'his article 

igure shows the geometry of the wind-tunnel model The 
configuration has a thin, low aspect ratio, highly swept wins 
mounted below the centerline of a slender body The wing if 
flat with a rounded leading edge. It should be noted that the 
exact wmgtip definition was not available and so the tip thick- 
ness was decreased to zero across three grid points 
Figure r shows the overview of the surface grid for the full- 

Te v rhefhne , ;°T 1 ! ,he f er,d ""f the *** ™ 

every other line). The reference length is taken from the mean 
erodynamic chord and the origin of the coordinates is set at 

SiSL M- ThC ,Wd I " CX,ended the d ~m 

boundary. The half-span grid used for the symmetric cases 


1 ln tn e streamwise direction. 1 16 points 

in the spanw.se direction, and 40 points normal to the body 
surface, for a total of 510,400 points. The bilateral symmetrJ 
ls 1 111 P os ed in the x-z plane at y = 0 (the center of 
e body). In the following computations, the grid is further 
divided into the upper and lower grids at the wing and the 

f-.ce OP Fo°r 8 fh eU f "n " ,s P rov,ded through a zonal inter- 
ace. For the full-span configuration used for the antisym- 
metric cases, the grid is mirrored to the other side (total of 

Tom ™ ,‘ h n the "T ber of grid points is dou bled 

Ini:,!:?!’ f Flow variables at the zonal interfaces were 
updated as soon as the adjoining zones were computed This 

fion! a , S r ,Mmpl,at Z T' in,erface for steady-state calcula- 

S , ini h T C Pr ° CCdarC Can be applied t0 unsteady cal- 
culations by alternating the sweep through zones at every time 

To treat the control surface movement without introducing 
additional zones, a small gap is introduced at the end of the 
control surface. This region is used to shear the grid when 

eZ , SU 1 f OSCillateS - The dynamic g fid ar °und a de- 
flected control surface was obtained by shearing every grid 

liTan” l ° C0ntr0 ' Wi * h the loca ' deflection. 

In the experiment, a transition strip was placed at the 15% 
chord. However, the report did not show significant differ- 
ences ,n comparisons of force measurement with and without 
he strip at the transonic regime.' In addition, the effect of 

de-o Thn" ,he t . Separation at «he leading edge was not very 
clear. Ihus, in the computation, a fully turbulent flow is as- 
sumcd. 

The grid lines on the body surface collapse to a point at 
the nose and extend upstream as a singular axis. The flow 
variables on the singular axis are given by taking an average 
from the surrounding grid points. When a computation starts 
impulsively from the freestream condition, the upwind method 
is not dissipative enough to damp the initial disturbances alo^ 
he axis. The central-difference option of the code was used 
o overcome th.s initial transient period. Since the upwind 
solution gave a crisper vortex structure for steady state, the 
upwind option was used for the rest of the calculations. 
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Results 

Steady Flap Deflection 

Figure 3 shows the steady pressures compare 
experiment at the 20. 50, and 80% semispanwise sections for 

f£Sff5 

S. to «* leading-edge »ooex The,, i » a 

minor discrepancy between the computation and the exper 

intent due to the difference in the loc ^ on ^ 
edf?e vortex The computation predicts the vortex at a siign y 
more inboard location than the experiment Possible sources 
of this difference are the effects of the transition strip and the 
wall of the wind tunnel. No data correction was applied to 
Either the computed or measured data. Overall, the computed 
result shows good agreement with the experiment. The pres- 
sure disti^^ on the body center also show good agree- 

m lt a Lt^n n dmg 8 result at the same flow condition with 
the flap deflected down by 8.3 deg “ g* Unwise 

Sn. although no «,uam»»c 

the trailing edge. The effect is not noticeable at the 20 A- 

section. The computed pressure profiles capture t e o* . 

tures well. The effect of the flap deflection is very small on 


Computation 

e Experiment, M»nro, •! »l. 





Fig. 3 Comparison of computed steady pressures with experiment; 
no flap deflection. 


Computation 

o Experiment, Manro, at at. 



Fig. 4 Comparison of computed steady pressures with experiment 
on the body surface; no flap deflection. 


Computation 

o Experiment, llanro, at al. 




x/c 


Fig. 5 Comparison of computed steady pressures with experiment, 
8.3-deg flap deflection. 


he body, and thus, the pressure distributions are not shown 

Table 1 shows the comparison of force coefficients. Both 
lormal force and moment coefficients show good ^ r «men 
vith the experimental data. The lift coefficients are 0.346 and 
) 310 with and without flap deflection, respectively. 

"Figure 6 shows the steady pressures for the full-span con- 
figuration with antisymmetric flap position. The right wing 
hfs the flap up by 8.3 deg and the left wing has the flap down 
hv 8 3 deg The pressure distributions show the largest di 
ference between The left and right wings at the 80% sect, on 
as expected. The plot also shows a discrepancy between the 
antisymmetric case and the symmetric case with the flap de- 
flected down. The antisymmetric position of the flaps ge 
erates a lower pressure above the left wing and a higher pres-^ 
sure above the right wing. This introduces a circulation around 
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Table I Comparison of force coefficients with experiment 


Computation 

Experiment' 


a = 7.93 deg. 

5 = 0 deg 

>y cy, 

0.298 -0.063 

0.295 - 0.065 


£* = 7.93 deg, 

S = 8.30 d eg 

Cmy 

0.332 -0.092 

0.328 -0.093 


Symmetric, flap down 
Anti-aym metric, flap down 
AntJ-symmetric, flap up 




C ? m P“ ris «n of computed steady pressures between the svm- 
metric and antisymmetric flap deflections. 


the x axis. Therefore, at the 80% section, the antisymmetric 
result shows smaller A C„ (and thus, smaller sectional lift) than 
the symmetric result. ; 

Figure 7 illustrates the structure of the vortical flowfield for 
the antisymmetric case by using the helicity density. Figure 
7a corresponds to the trailing edge. Note that the helicity 
density was computed on the grid points, and thus, the cross- 

not i exactiy the p ,ane at a c ° nsta nt * 

location. The other plots are taken downstream at intervals 
of approximately 0.03 and 0.06c. The primary and secondary 

can bareN n h f H d ° VCr b °‘ h wings ' The win g‘'P vortices 
can barely be found next to the secondary vortices. An in- 

terestmg feature is that the wake vortex sheet shed from the 

As P fh7f| 0n a UP !T d merges with the secondary vortex. 

do ser to 6 Z ^ ^ ^’“P V ° rteX sheet Torres 

closer to the primary vortex. Due to the interaction by the 

^ '.' '7 vortex - th e rolled-up vortex is displaced towards the 
-ah ndary vo [ ,ex - and thus, they merge more quickly. On the 
other hand, the vortex interaction is moderate when the flap 
is down. Instead, the vortex itself is stronger because of the 
camber introduced by the deflected flap. Without the flap 
deflection the corresponding wake structure basically falls 
into the middle of the left and right wake structures ’ 



uens,, y contour plots in the wake for the 
down”™™ CaSe: a> ' railin8 M,8e ’ b ’ ° °' V downs ‘>-cam. and cl 0.06c 



P: Primary vortex 
S: Saoondary vortex 
R: Rolted-up wake 
F: Rap-lip vortex 



The density contour plot in the crossflow plane (the true 
.W plane) at x = 2.6 is shown in Fig. 8. The right half of the 

the left r haT ndS t0 the . wake f ° r the u P ward tla P position and 
the eft half corresponds to the wake for the downward flap 

position. On the left-hand side (LHS), four low-density re P 

vone'xT 3 b | e id d; 1} “T pnmary vortex p ' 2) secondar y 
vortex 5,. 3) rolled-up wake vortex sheet R, and 4) the flap 

inboard-tip vortex F. The wingtip vortex is weaker than these 

tex rt and rm dearlyobserved here - The secondary vor- 
tex and the rolled-up wake vortex sheet are really the same 

vorbcal region and they merge rapidly as shown in Fig. 7 
The flap inboard-tip vortex can be seen only on the left side 
“ fla r ls defIected down On the right-hand side 

height of H ^ , OW denS " y reg, ° ns are found - Comparing the 
height of the two primary vortices, the left one is located 

lower due to the flap deflected down. Also, the lower density 
at the center of the left primary vortex indicates a stronger 
■ ICI y ' Wben we dr aw a line connecting the center of the 
pnmary vortex and the center of the rolled-up wake vortex 

LHs’whd rh makCS 3 3? ‘ deg a " gle *° the vertical on the 
the RHS if 6 c ^ res Po nd >ng hne makes a 44-deg angle on 

Hill H h Th hC rolled ' u P vortex sheet on the RHS is 
displaced more toward the secondary vortex. 

Oscillatory Flap Motion 

The capability of ENSAERO code used in this work to 
compute unsteady flows with oscillating flaps is previously 
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have slightly higher lift on the average (the lift responses are 
shown later). 

Unsteady pressure responses to the flap motion are shown 
in magnitude and phase angle in Fig. 10. In the magnitude 
plots, the left and right wings for the antisymmetric case give 
identical responses again because the only difference in the 
flap motion is the 180-deg phase angle. Also, antisymmetric 
results show slightly lower responses than the symmetric result 
because the 180-deg phase difference cancels the pressure 
variation at the body centerline. At the 80% spanwise section, 
a sharp peak is found at the hinge line of the control surface, 
which is located at the 75% chord. After the 90% chord, the 
magnitude of C p drops and rises again near the trailing edge. 
A similar trend can be found at the 50% spanwise section. 
In the phase plots, the left and right wings for the antisym- 
metric case clearly show a 180-deg phase difference at the 
80% section. Corresponding to the magnitude plot at this 
section, a jump of the phase change by 180 deg can be found 
around the 90% chord. This unsteady response is related to 
the interaction between the primary vortex and the wake 
vortex sheet as explained later. 

At the inboard sections, the phase plots for the antisym- 
metric case are not as smooth as those for the symmetric case. 
The 180-deg phase difference between the left and right wings 
is not seen, either. Since the magnitude for the antisymmetric 
case is nearly zero in this region, the phase plots for the 
antisymmetric case are more sensitive to numerical errors. 
Second-order time accuracy may be needed to achieve a per- 
fect 180-deg phase difference for reasonable time step sizes. 
The result was not improved very much by simply halving the 
time step size. For practical purposes, the current result is 
accurate enough because the instantaneous C p plots on the 
left wing still coincide with those of the right wing after a half 
cycle (180-deg difference). 

Figure 11 shows the instantaneous density contour plots in 
the crossflow plane (true y-z plane) at x = 2.6, similar to Fig. 

8. Five instances were chosen at 0-, 45-, 90-, 135-, and 180- 
deg phase angles. The left wing starts from 0 deflection, goes 
down, and comes back to 0 deflection. The right flap starts 
from 0 deflection, goes up, and comes back to 0 deflection. 
The rest of the half cycle is a mirror image of this half cycle. 
By comparing Figs. 1 la and 1 lc, the left primary vortex goes 
down and becomes stronger as the flap goes down and the 
right primary vortex goes up and becomes weaker as the flap 
goes up. 

When we draw lines connecting the center of the primary 
vortex and the center of the rolled-up wake vortex sheet, the 
angles that the line on the RHS make with the vertical line 
vary more than the corresponding angles on the LHS. This 
indicates that the primary vortex has a stronger influence on 
the wake vortex sheet from the flap when it goes up. This 
explains the pressure response near the trailing edge at the 
80% section observed earlier in Fig. 10. It also explains the 
180-deg phase jump because the disturbance occurs when the 
flap is up, which is negative deflection-by-definition. 

By comparing Figs. 1 lb and 1 Id, a slight discrepancy can 
be found in the vortex structure. This suggests hysteresis in 
the lift response. Figure 12 shows the lift responses with re- 
spect to the flap deflection angle of the right wing. As ob- 
served in Fig. 11, the right wing shows the hysteresis as the 
flap is deflected up (negative deflection). The total lift re- 
sponse shows that the unsteady case has more lift than the 
steady case {€, = 0.310). 

The symmetric case shows the hysteresis when the flap is 
deflected down. This increase in the lift is associated with the 
primary vortex enhanced by the downward deflection of the 
flap. The overall lift for the symmetric case is higher than the 
antisymmetric case. This is consistent to the steady pressure 
result shown in Fig. 6. On the other hand, the antisymmetric 
case introduces a circulation around the x axis as discussed 
earlier in the steady results. This circulation reduces the strength 


P: Primary vortex 
S: Secondary vortex 
R: Rolled-up wake 
F: Flap-tip vortex 



x = 2.6 





Fig. 11 Crossflow density contour plots at x = 2.6 for the antisym- 
metric oscillatory motion of the flaps. 
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Right flap daflaction (dag) 

Fig. 12 Comparison of lift responses between the symmetric and 
antisymmetric oscillatory motions of the flaps. 

of the primary vortex with the flap deflected down (the stronger 
side) but increases the strength with the flap up (the weaker 
side)’ As a result, the hysteresis appears during the flap de- 
flected up in the antisymmetric case. 

Conclusions 

Unsteady Navier-Stokes simulations of vortical flows in 
the transonic regime over a rigid arrow-wing configuration 
with trailing-edge control surfaces have been performed. 
Computations have been made at a moderate angle of attack, 
where a leading-edge separation occurs. Computed steady 
pressures and integrated force coefficients with and without 


control surface deflection show good agreement with the ex- 
perimental data. 

A comparison of response characteristics between the sym- 
metric and antisymmetric control surface oscillatory motions 
on the left and right wings is also presented. For unsteady 
computations, the grid moves every time step following the 
deflection of the control surface. The antisymmetric case pro- 
duces higher lift than the steady case with no deflection, and 
the symmetric case produces higher lift than the antisymmetric 
case. The unsteady pressure responses indicate strong influ- 
ence of the flap motion at the flap trailing edge for both 
symmetric and antisymmetric cases. The detailed analysis of 
the wake structure revealed a strong interaction between the 
primary vortex and the wake vortex sheet from the flap re- 
gion. 
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